Search for: All records

Abstract Spatial patterns of stellar elemental abundances encode rich information about a galaxy’s formation history. We analyze the radial, vertical, and azimuthal variations of metals in stars, both today and at formation, in the FIRE-2 cosmological simulations of Milky Way (MW)-mass galaxies, and we compare them with the MW. The radial gradient today is steeper (more negative) for younger stars, which agrees with the MW, although radial gradients are shallower in FIRE-2. Importantly, this age dependence was present already at birth: radial gradients today are only modestly (≲0.01 dex kpc⁻¹) shallower than at birth. Disk vertical settling gives rise to negative vertical gradients across all stars, but vertical gradients of mono-age stellar populations are weak. Similar to the MW, vertical gradients in FIRE-2 are shallower at larger radii, but they are overall shallower in FIRE-2. This vertical dependence was present already at birth: vertical gradients today are only modestly (≲0.1 dex kpc⁻¹) shallower than at birth. Azimuthal scatter is nearly constant with radius, and it is nearly constant with age ≲8 Gyr ago but increases for older stars. Azimuthal scatter is slightly larger (≲0.04 dex) today than at formation. Galaxies with larger azimuthal scatter have a stronger radial gradient, implying that azimuthal scatter today arises primarily from the radial redistribution of gas and stars. Overall, spatial variations of stellar metallicities show only modest differences between formation and today; spatial variations today primarily reflect the conditions of stars at birth, with spatial redistribution of stars after birth contributing secondarily. 

ABSTRACT We analyse the stellar distributions on the [Fe/H]–[Mg/Fe] plane for 11 Milky Way-mass galaxies from the FIRE-2 cosmological baryonic zoom-in simulations. Alpha-element bimodality, in the form of two separate sequences on the [Fe/H]–[Mg/Fe] plane, is not a universal feature of disc galaxies. Five galaxies demonstrate double sequences with the $$\alpha$$-enriched one being older and kinematically hotter, in qualitative agreement with the high-$$\alpha$$ and low-$$\alpha$$ populations in the Milky Way disc; three galaxies have unimodal distribution, two show weakly bimodal features where low-$$\alpha$$ sequence is visible only over a short range of metallicities, and one show strong bimodality with a different slope of high-$$\alpha$$ population. We examine the galaxies’ gas accretion history over the last 8 Gyr, when bimodal sequences emerge, and demonstrate that the presence of the low-$$\alpha$$ sequence in the bimodal galaxies is related to the recent infall of metal-poor gas from the circumgalactic medium that joins the galaxy in the outskirts and induces significant growth of the gas discs compared to their non-bimodal counterparts. We also analyse the sources of the accreted gas and illustrate that both gas-rich mergers and smooth accretion of ambient gas can be the source of the accreted gas, and create slightly different bimodal patterns. 

The stellar halos of galaxies, primarily formed through the accretion and merger of smaller objects, are an important tool for understanding the hierarchical mass assembly of galaxies. However, the inner regions of stellar halos in disk galaxies are predicted to have an in situ component that is expected to be prominent along the major axis. Kinematic information is crucial to disentangle the contribution of the in situ component from the accreted stellar halos. The low surface brightness of stellar halos makes it inaccessible with traditional integrated light spectroscopy. In this work, we used a novel technique to study the kinematics of the stellar halo of the edge-on galaxy NGC 4945. We couple new deep Multi Unit Spectroscopic Explorer spectroscopic observations with existingHubbleSpace Telescope imaging data to spectroscopically measure the line-of-sight (LOS) heliocentric velocity and velocity dispersion in two fields at a galactocentric distance of 12.2 kpc (outer disk field) and 34.6 kpc (stellar halo field) along the NGC 4945 major axis, by stacking individual spectra of red giant branch and asymptotic giant branch stars. We obtained a LOS velocity and dispersion of 673 ± 11 km s⁻¹and 73 ± 14 km s⁻¹, respectively, for the outer disk field. This is consistent with the mean HI velocity of the disk at that distance. For the halo field, we obtained a LOS velocity and dispersion of 519 ± 12 km s⁻¹and 42 ± 22 km s⁻¹. The halo fields’ velocity measurement is within ∼40 km s⁻¹from the systemic LOS velocity of NGC 4945, which is 563 km s⁻¹, suggesting that its stellar halo at 34.6 kpc along the major axis is counter-rotating and its origins are likely to be the result of accretion. This provides the first-ever kinematic measurement of the stellar halo of a Milky Way-mass galaxy outside the Local Group from its resolved stellar population. Thus, we have established a powerful technique for measuring the velocity field for the stellar halos of nearby galaxies. 

Abstract Mergers of and interactions between galaxies imprint a wide diversity of morphological, dynamical, and chemical characteristics in stellar halos and tidal streams. Measuring these characteristics elucidates aspects of the progenitors of the galaxies we observe today. The M81 group is the perfect galaxy group to understand the past, present, and future of a group of galaxies in the process of merging. Here, we measure the end of star formation (t₉₀) and metallicity ([M/H]) of the stellar halo of M82 and the eastern tidal stream of NGC 3077 to: (1) test the idea that M82 possesses a genuine stellar halo, formed before any interaction with M81; (2) determine if NGC 3077's tidal disruption is related to the star formation history in its tails; and (3) create a timeline of the assembly history of the central trio in the M81 group. We argue that M82 possesses a genuine, metal-poor ([M/H] ∼ −1.62 dex) stellar halo, formed from the merger of a small satellite galaxy roughly 6.6 Gyr ago. We also find that the stars present in NGC 3077's tails formed before tidal disruption with M81, and possess a roughly uniform metallicity as shown in S. Okamoto et al., implying that NGC 3077's progenitor had significant population gradients. Finally, we present a timeline of the central trio’s merger/interaction history. 

Abstract In the Λ-Cold Dark Matter model of the universe, galaxies form in part through accreting satellite systems. Previous works have built an understanding of the signatures of these processes contained within galactic stellar halos. This work revisits that picture using seven Milky Way–like galaxies in the Latte suite of FIRE-2 cosmological simulations. The resolution of these simulations allows a comparison of contributions from satellites above M * ≳ 10 × 7 M ⊙ , enabling the analysis of observable properties for disrupted satellites in a fully self-consistent and cosmological context. Our results show that the time of accretion and the stellar mass of an accreted satellite are fundamental parameters that in partnership dictate the resulting spatial distribution, orbital energy, and [ α /Fe]-[Fe/H] compositions of the stellar debris of such mergers at present day. These parameters also govern the resulting dynamical state of an accreted galaxy at z = 0, leading to the expectation that the inner regions of the stellar halo ( R GC ≲ 30 kpc) should contain fully phase-mixed debris from both lower- and higher-mass satellites. In addition, we find that a significant fraction of the lower-mass satellites accreted at early times deposit debris in the outer halo ( R GC > 50 kpc) that are not fully phased-mixed, indicating that they could be identified in kinematic surveys. Our results suggest that, as future surveys become increasingly able to map the outer halo of our Galaxy, they may reveal the remnants of long-dead dwarf galaxies whose counterparts are too faint to be seen in situ in higher-redshift surveys. 

ABSTRACT Milky Way-mass galaxies in the FIRE-2 simulations demonstrate two main modes of star formation. At high redshifts star formation occurs in a series of short and intense bursts, while at low redshifts star formation proceeds at a steady rate with a transition from one mode to another at times ranging from 3 to 7 Gyr ago for different galaxies. We analyse how the mode of star formation affects iron and alpha-element abundance. We find that the early bursty regime imprints a measurable pattern in stellar elemental abundances in the form of a ‘sideways chevron’ shape on the [Fe/H] – [O/Fe] plane and the scatter in [O/Fe] at a given stellar age is higher than when a galaxy is in the steady regime. That suggests that the evolution of [O/Fe] scatter with age provides an estimate of the end of the bursty phase. We investigate the feasibility of observing of this effect by adding mock observational errors to a simulated stellar survey and find that the transition between the bursty and steady phase should be detectable in the Milky Way, although larger observational uncertainties make the transition shallower. We apply our method to observations of the Milky Way from the Second APOKASC Catalogue and estimate that the transition to steady star formation in the Milky Way happened 7 – 8 Gyrs ago, earlier than transition times measured in the simulations. 

Abstract The shape and orientation of dark matter (DM) halos are sensitive to the microphysics of the DM particles, yet in many mass models, the symmetry axes of the Milky Way’s DM halo are often assumed to be aligned with the symmetry axes of the stellar disk. This is well motivated for the inner DM halo, but not for the outer halo. We use zoomed-in cosmological baryonic simulations from the Latte suite of FIRE-2 Milky Way–mass galaxies to explore the evolution of the DM halo’s orientation with radius and time, with or without a major merger with a Large Magellanic Cloud analog, and when varying the DM model. In three of the four cold DM halos we examine, the orientation of the halo minor axis diverges from the stellar disk vector by more than 20° beyond about 30 galactocentric kpc, reaching a maximum of 30°–90°, depending on the individual halo’s formation history. In identical simulations using a model of self-interacting DM withσ= 1 cm²g⁻¹, the halo remains aligned with the stellar disk out to ∼200–400 kpc. Interactions with massive satellites (M≳ 4 × 10¹⁰M_⊙at pericenter;M≳ 3.3 × 10¹⁰M_⊙at infall) affect the orientation of the halo significantly, aligning the halo’s major axis with the satellite galaxy from the disk to the virial radius. The relative orientation of the halo and disk beyond 30 kpc is a potential diagnostic of self-interacting DM, if the effects of massive satellites can be accounted for. 

ABSTRACT We characterize the 3D spatial variations of [Fe/H], [Mg/H], and [Mg/Fe] in stars at the time of their formation, across 11 simulated Milky Way (MW)- and M31-mass galaxies in the FIRE-2 simulations, to inform initial conditions for chemical tagging. The overall scatter in [Fe/H] within a galaxy decreased with time until $$\approx 7 \, \rm {Gyr}$$ ago, after which it increased to today: this arises from a competition between a reduction of azimuthal scatter and a steepening of the radial gradient in abundance over time. The radial gradient is generally negative, and it steepened over time from an initially flat gradient $$\gtrsim 12 \, \rm {Gyr}$$ ago. The strength of the present-day abundance gradient does not correlate with when the disc ‘settled’; instead, it best correlates with the radial velocity dispersion within the galaxy. The strength of azimuthal variation is nearly independent of radius, and the 360 deg scatter decreased over time, from $$\lesssim 0.17 \, \rm {dex}$$ at $$t_{\rm lb} = 11.6 \, \rm {Gyr}$$ to $$\sim 0.04 \, \rm {dex}$$ at present-day. Consequently, stars at $$t_{\rm lb} \gtrsim 8 \, \rm {Gyr}$$ formed in a disc with primarily azimuthal scatter in abundances. All stars formed in a vertically homogeneous disc, Δ[Fe/H]$$\le 0.02 \, \rm {dex}$$ within $$1 \, \rm {kpc}$$ of the galactic mid-plane, with the exception of the young stars in the inner $$\approx 4 \, \rm {kpc}$$ at z ∼ 0. These results generally agree with our previous analysis of gas-phase elemental abundances, which reinforces the importance of cosmological disc evolution and azimuthal scatter in the context of stellar chemical tagging. We provide analytic fits to our results for use in chemical-tagging analyses. 

Abstract It is not yet settled how the combination of secular processes and merging gives rise to the bulges and pseudobulges of galaxies. The nearby (D∼ 4.2 Mpc) disk galaxy M94 (NGC 4736) has the largest pseudobulge in the local universe, and offers a unique opportunity for investigating the role of merging in the formation of its pseudobulge. We present a first ever look at M94's stellar halo, which we expect to contain a fossil record of M94's past mergers. Using Subaru's Hyper Suprime-Cam, we resolve and identify red giant branch (RGB) stars in M94's halo, finding two distinct populations. After correcting for completeness through artificial star tests, we can measure the radial profile of each RGB population. The metal-rich RGB stars show an unbroken exponential profile to a radius of 30 kpc that is a clear continuation of M94's outer disk. M94's metal-poor stellar halo is detectable over a wider area and clearly separates from its metal-rich disk. By integrating the halo density profile, we infer a total accreted stellar mass of ∼2.8 × 10⁸M_⊙, with a median metallicity of [M/H] = −1.4. This indicates that M94's most-massive past merger was with a galaxy similar to, or less massive than, the Small Magellanic Cloud. Few nearby galaxies have had such a low-mass dominant merger; therefore we suggest that M94's pseudobulge was not significantly impacted by merging.